Observation of Phase Fluctuations in Bose-Einstein Condensates 
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The occurrence of phase fluctuations due to thermal excitations in Bose-Einstein condensates (BECs) 
is studied for a variety of temperatures and trap geometries. We observe the statistical nature 
of the appearence of phase fluctuations and characterize the dependence of their average value on 
temperature, number of particles and the trapping potential. We find pronounced phase fluctuations 
for condensates in very elongated traps in a broad temperature range. The results are of great 
importance for the realization of BEC in quasi ID geometries, for matter wave interferometry with 
BECs, as well as for coherence properties of guided atom laser beams. 
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Since the first experimental realization of Bose- 
Einstein condensation in dilute atomic gases (lj, there 
has been enormous interest in the coherence properties 
of BECs. In particular, the phase coherence is essential 
for applications of BEC in matter wave interferometry; 
it also sets limits on the coherence of atom lasers, and 
guided atom laser beams. For a trapped 3D condensate 
well below the BEC transition temperature T c , recent 
experiments have confirmed the phase coherence, e.g., it 
was shown that the coherence length is equal to the con- 
densate size [||J|]. However, phase coherence is not an 
obvious property of BEC. In particular, it is expected 
that low-dimensional (ID and 2D) quantum gases dif- 
fer qualitatively from the 3D case in this respect || 
Recently, it was shown theoretically § that for very elon- 
gated condensates phase fluctuations can be pronounced 
already in the equilibrium state of the usual 3D ensem- 
ble, where the density fluctuations are suppressed. The 
phase coherence length in this case can be smaller than 
the axial size of the sample. This is referred to as the 
regime of quasicondensation ||. The detailed charac- 
terization of phase fluctuations in condensates is thus of 
great importance for applications of BEC, especially for 
recent attempts to reach BEC in elongated micro-circuit 
geometries . The temperature dependence of the co- 
herence of an atom laser beam was studied in fl3| ] . 

In this Letter we report on systematic studies of BEC 
of 87 Rb atoms in the regime of phase-fluctuating con- 
densates. We achieved this regime in highly anisotropic 
traps leading to a strongly elongated shape of the con- 
densate. We observe the phase fluctuations by measuring 
the density distribution of the released cloud after bal- 



listic expansion. By varying the temperature and the 
aspect ratio of the trapping potential, we study the con- 
tinuous transition from the usual 3D regime, where phase 
fluctuations of the condensate are low, into the regime 
of strong phase fluctuations. We show that a phase co- 
herent matter wave is not a direct outcome of BEC in 
strongly elongated geometries but can only be achieved 
for very low temperatures, well below the BEC transition 
temperature T c . 

Fluctuations of the phase of a Bose condensate are re- 
lated to thermal excitations and always appear at finite 
temperature. However, as shown in Ref. the fluctu- 
ations depend not only on temperature but also on the 
trap geometry, and on the particle number. Typically, 
fluctuations in spherical traps are strongly suppressed as 
the wavelengths of excitations are smaller than the size 
of the atomic cloud. In contrast, wavelengths of the exci- 
tations in strongly elongated traps can be larger than the 
radial size of the cloud. In this case, the low-energy axial 
excitations acquire a ID character and can lead to more 
pronounced phase fluctuations, although the density fluc- 
tuations of the equilibrium state are still suppressed. Due 
to the pronounced phase fluctuations the coherence prop- 
erties of elongated condensates can be significantly al- 
tered as compared with previous observations. In par- 
ticular, the axial coherence length can be much smaller 
than the size of the condensate, which can have dramatic 
consequences for practical applications. 

The experiment was performed with Bose-Einstein 
condensates of up to iVo = 5 x 10 5 rubidium atoms in 
the \F = 2,rriF = +2) state of a cloverleaf-type magnetic 
trap. Further details of our apparatus were described 
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previously |14||n|. The fundamental frequencies of the 
magnetic trap are lo x = 2ir x 14 Hz and lo p = 2tt x 365 Hz 
along the axial and radial direction, respectively. Due 
to the highly anisotropic confining potential with an as- 
pect ratio A = ui p /w x of 26, the condensate is already 
elongated along the horizontal x axis. In addition, fur- 
ther radial compression of the ensemble by means of a 
superimposed blue detuned optical dipole trap is possi- 
ble |16|. We performed the measurements for radial trap 
frequencies ui p between 2ir x 138 Hz and 2ir x 715 Hz cor- 
responding to aspect ratios A between 10 and 51. After rf 
evaporative cooling to the desired temperature, we wait 
for 1 sec (with rf 'shielding') to allow the system to reach 
an equilibrium state. We then switch off the trapping po- 
tential within 200 [is and allow a variable time-of- flight. 



velocity fields provided by the phase fluctuations into 
modulations of the density. However, after switching 
off the trap, the mean-field interaction rapidly decreases 
and the axial velocity fields are then converted into the 
density distribution. We have performed numerical sim- 
ulations of the 3D Gross-Pitaevskii equation to under- 
stand quantitatively how phase fluctuations lead to the 
build up of stripes in the density distribution. We as- 
sume that initially the condensate had an equilibrium 
density profile, and a random fluctuating phase <fi(x). 
For elongated BECs @] the phase can be represented as 
4>(x) = Ejli <l>j( x ), where 
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FIG. 1. Absorption images and corresponding density pro- 
files of BECs after 25 ms time-of-flight taken for various tem- 
peratures T and aspect ratios [A = 10 (a), 26 (b), 51 (c)]. 

Figure ^ shows examples of experimental data for var- 
ious temperatures T < T c and aspect ratios A. The usual 
anisotropic expansion of the condensate related to the 
anisotropy of the confining potential is clearly visible in 
the absorption images. The line density profiles below 
reflect the parabolic shape of the BEC density distribu- 
tion. Remarkably, we observe pronounced stripes in the 
density distribution in some cases shown in Fig. ^. On 
average these stripes are more pronounced at high aspect 
ratios of the trapping potential [Fig. 0(c)], high temper- 
atures (bottom row of Fig. ^), and low atom numbers. 

The appearence of stripes can be understood qualita- 
tively as follows. As mentioned above, within the equi- 
librium state of a BEC in a magnetic trap the density 
distribution remains largely unaffected even if the phase 
fluctuates fa]. The reason is that the mean- field inter- 
particle interaction prevents the transformation of local 



Here 6j = hu} x + 3)/4 is the spectrum of low-energy 
axial excitations |L7|], Pj 1 are Jacobi polynomials, g = 
Airh 2 a/m, a > is the scattering length, and R (L) is the 
radial (axial) condensate size. In Eq. ([!]), quasiparticle 
creation and annihilation operators have been replaced 



by complex amplitudes ctj and 



To reproduce the 



quantum statistical properties of the phase, ay and a* 
were sampled as random variables with a zero mean value 



and (\oej 



Nj , where Nj 



[exp(ej/k B T) - 1] _1 is the 



occupation number for the quasiparticle mode j. 
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FIG. 2. (a) A typical initial phase distribution calculated 
for w x = 2tt x 14Hz, u p = 2tt x 508 Hz, JV = 2x 10 5 , and 
T — 0.5 T c . (b) Corresponding density profile after 25 ms 
time-of-flight from simulations (solid line) and analytical the- 
ory (dashed line). 

Alternatively, the formation of the stripes has been 
studied by using the local density approximation for the 
axial profile and by relying on the scaling approach |lq ] 
for the radial expansion of the cloud. The rescaled Gross- 
Pitaevskii equation was then linearized with respect to 
small fluctuations of the phase gradients and the den- 
sity. Our method accounts for the transformation of ini- 
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tially phonon-like excitations into particle-like ones and 
covers both the (initially) hydrodynamic and the (ulti- 
mately) free-particle regimes of expansion. For a con- 
densate which initially is in the Thomas-Fermi regime 
and has a chemical potential /i, the radially integrated 
density fluctuations Sn(x) after a time-of- flight t in the 
interval h/Tilu 2 , t 3> [ijfvJ 1 are given by 
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Here the profile n${x) is the radially integrated density 
for the unperturbed condensate. From Eq. (j^) one ob- 
tains a closed relation for the mean square density fluc- 
tuations (<7 B ec/«-o) 2 by averaging (Sn/no) 2 over different 
realizations of the initial phase. In the central part of the 
cloud (x « 0) we find 
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where r = u p t. In Eq. (g) we introduced a characteristic 
temperature k B T^, = l5(hui x ) 2 N /32fi. For T$ < T c one 
expects the regime of quasicondensation for the initial 
cloud in the temperature interval < T < T Cl i.e., 
the regime where the coherence length is smaller than 
the condensate size [|| . The analytical expressions agree 
very well with the numerical simulations (see Fig. ||). 
Note that Eq. (||) provides a direct relation between the 
observed density fluctuations and the temperature, and 
thus can be used for thermometry at very low T. 
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FIG. 3. Measurement of Obec/io versus temperature in 
four different trap geometries. The dotted lines represent the 
average detection noise or /no- All data was taken for T < T c . 



To determine the amount of phase fluctuations exper- 
imentally, we systematically studied the formation and 
structure of stripes in the atomic density distribution as 
a function of the trapping potential and temperature. 
For each realization of a BEC, the observed density dis- 
tribution was integrated along the radial direction and 
then fitted by a bimodal function with the integrated 
parabolic Thomas-Fermi distribution for the condensate 
fraction and a Gaussian for the thermal cloud. For each 
image we obtained standard deviations ct B ec of the ex- 
perimental data from the fit in the central region of the 
condensate fraction (half width of full size). 

The temperature dependence of a BEC is shown in Fig. || 
for four different trap configurations with aspect ratios of 
A = 10, 26, 36 and 51. To account for shot-to-shot varia- 
tions in the atom number, the standard deviations were 
normalized to the fitted peak density uq in the conden- 
sate. The dotted lines show the average detection noise 
crr/no for data from each trap, obtained from the ther- 
mal wings of the cloud. 

Since the initial phase of a Bose condensate is mapped 
into its density distribution after time-of- flight, the quan- 
tity <7 BE c is a direct measure of the initial phase fluc- 
tuations. Note, however, that this method reflects the 
instantaneous phase of the BEC at the time of release 
and, therefore, images taken at the same initial condi- 
tions can look significally different. Indeed, we observe a 
large spread of our experimental data (see Fig. |^), which 
clearly demonstrates the statistical character of the phase 
fluctuations. Furthermore, our data characterizes phase 
fluctuations for various trap geometries. For A — 10 the 
data falls into our detection noise and, therefore, the 
phase fluctuations are hardly observable. Considerably 
different images are observed for more elongated traps 
with A = 26,36 and 51. In these cases we observe sig- 
nificant deviations from the Thomas-Fermi distribution 
at high temperatures, which is a clear indication of the 
presence of phase fluctuations. For lower temperatures, 
reached by further evaporation, phase fluctuations be- 
come reduced due to both, the reduced excitation spec- 
trum at lower temperature and the increasing number of 
atoms in the condensate fraction. However, for even fur- 
ther evaporation the number of Bose-condensed atoms 
iVo decreases leading to phase fluctuations in agreement 
with Eq. (H). In fact, it was very difficult to reduce the 
density modulation below the noise limit by rf evapora- 
tion in the case of our tightest trap with A = 51. 

As observed in our experiment phase fluctuations for 
BECs in elongated traps continuously decrease with re- 
ducing T and increasing iVo- It is not possible to de- 
termine a cut-off for the phase fluctuations, rather they 
decrease until they cannot be resolved below our noise 
limit. Hence all experiments with BECs at finite tem- 
perature in tightly confining elongated potentials will be 
subject to axial phase fluctuations. 

In order to obtain general information about the phase 
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fluctuations, we average the observed standard deviations 
over many realizations within a small temperature and 
particle number interval. This allows us to use Eq. (||) 
and to compare directly the averaged measured value of 
[K E c/no) 2 ]!/ p 2 to the predicted value of (o- BEC /n )thcor y 
for various values of T, Nq,lu x ,lu p , and t (see Fig. ^). 
The theoretical value takes also the limited experimental 
imaging resolution into account. 
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FIG. 4. Average standard deviation of the measured line 
densities [ (a B Ec/no) 2 ]Hp compared to the theoretical value 
of (cr B Ec/wo)tht!ory obtained from Eq. (^). The dashed line is 
a lit to the experimental data. 

Our experimental results follow the expected general 
dependence very well. With the direct link of the phase 
fluctuations in the magnetic trap to the observed density 
modulation given by Eq. (|^), our data therfore confirms 
the predicted general behavior of phase fluctuations in 
elongated BECs. However, the measured values are ap- 
proximately by a factor of 2 smaller than those predicted 
by theory. This discrepancy could be due to, e.g., a re- 
duction of the observed contrast caused by a small tilt 
in the detection laser beam with respect to the radial 
stripes. Most of the observed experimental data, which 
exhibit fluctuations well above the noise limit, correspond 
to T > T^. This implies that the measurements were 
performed in the regime of quasicondensation, i.e., the 
phase coherence length 1$ = LT^/T of the initial con- 
densate was smaller than the axial size L. For instance, 
for A = 51, T = 0.5 T c , and N = 3 x 10 4 , one obtains 
h « L/3. 

In conclusion, we have presented detailed experimental 
and theoretical studies of a BEC state with fluctuating 
phase. The flexible method of ballistic expansion allows 
to measure phase fluctuations under various experimental 
conditions, especially for various trap geometries. Bragg 
spectroscopy can also be used to measure the influence of 
phase fluctuations on the momentum distribution of par- 
ticles in the axial direction, but this method is difficult 
to apply to very elongated traps. In addition, the ballis- 
tic expansion for very long times can visualize extremely 
small phase fluctuations not accessible to the usual res- 
olution of Bragg spectroscopy. By measuring the phase 



fluctuations and comparing the temperature with we 
have demonstrated instances, where the phase coherence 
length was smaller than the axial size of the condensate, 
i.e., the initial cloud was in the quasicondensate state. 
Our results set severe limitations on applications of BECs 
in interferomctric measurements, and for guided atom 
laser beams. Our experimental method combined with 
the theoretical analysis [Eq. (0)] provides a method of 
BEC thermometry. Further studies of quasicondensates, 
e.g., with respect to superfluid properties, will be neces- 
sary to obtain a full understanding of the phase coherence 
properties of ultracold atomic gases. 
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